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ABSTRACT 

We show that several features reminiscent of short-hard Gamma-ray Bursts (GRBs) arise naturally 
when Quark-Novae occur in low-mass X-ray binaries born with massive neutron stars (> I.GMq) 
and harboring a circumbinary disk. Near the end of the first accretion phase, conditions are just 
right for the explosive conversion of the neutron star to a quark star (Quark- Nova). In our model, 
the subsequent interaction of material from the neutron star's ejected crust with the circumbinary 
disk explains the duration, variability and near-universal nature of the prompt emission in short- 
hard GRBs. We also describe a statistical approach to ejecta break-up and collision to obtain the 
photon spectrum in our model, which turns out remarkably similar to the empirical Band function 
(|Band et al.lll993D . We apply the model to the fluence and spectrum of GRB 000727, GRB 000218, 
and GRB980706A obtaining excellent fits. Extended emission (spectrum and duration) is explained 
by shock-heating and ablation of the white dwarf by the highly energetic ejecta. Depending on the 
orbital separation when the Quark-Nova occurs, we isolate interesting regimes within our model when 
both prompt and extended emission can occur. We find that the spectrum can carry signatures typical 
of Type Ib/c SNe although these should appear less luminous than normal type Ib/c SNe. Late X- 
ray activity is due to accretion onto the quark star as well as its spin-down luminosity. Afterglow 
activity arise from the expanding shell of material from the shock-heated expanding circumbinary 
disk. We find a correlation between the duration and spectrum of short-hard GRBs as well as modest 
hard-to-soft time evolution of the peak energy. 

Subject headings: Stars: evolution, stars: binary, stars: neutron, supernovae: general, gamma-ray 
burst: general 



1. INTRODUCTION 

Intensive sampling of gamma-ray burst (GRB) light 
curves by Swift has revealed that at least in some short 
GRBs (henceforth shGRBs), prompt emission is fol- 
lowed by softer, extended emission lasting tens to hun- 
dreds of seco nds. For example , the lightcurve of GRB 
050709 ( Villasenor et all I2005D has a short-hard pulse 
Tsn ~ 0.2s and a long-so ft pulse Tgo - 130s; GRB 050724 
(jBarthelmv et al.l 120051 ) has prominent emission lasting 
for ~ 3 s followed by a long, soft, less prominent emis- 
sion peaking at ~ 100 s after the trigger, while XRT 
observations reveal strong flare-like activities within the 
first hundreds of seconds. In GRB 060614, the lightcurve 
has a short-hard episode followed by an extended soft 
emission component with strong spectral evolution. On 
average, the principal properties of the intense prompt 
component (< 2 s) of short bursts with and without 
extended emission (EE) are indisting uishable, suggest- 
ing t hat short bursts are "universal" (|Norris fc Gehrelsl 
[20081) . While the origin of such extended emission (seen 
roughly in a quarter of Swift shGRBs) is debated, it can 
provide some clues to the identity of the elusive mecha- 
nism of shGRBs. 

The association of some shGRBs with early type galax- 
ies that have low star formation rates suggests that 
short bursts arise from a different progenitor mecha- 
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nism than long bursts and that neutron star-black hole 
(NS-BH) or double neutron star (NS-N S) mergers might 
be at play (iBloomfcProc hasfe] 120061 ). But how does 
one obtain extended emission in this case? Extended 
accreti on episodes in NS-BH mergers have been pro - 
posed (|Barthelmv et al.l I2005D . while iRosswod (|2006f ) 
suggest that some debris may be launched during the 
merger process, which would fall back later to power 
flares at late time s. Alternatively, disk fragmentation 
(jPerna et al.l l2006f) or magneti c field barrier near the 
accretor ( Proga fc Zhang] |2006[ ) may induce intermit- 
tent accretion that power the flares. Other types of 
mergers such a s a white dwarf-neutron star (WD-NS) 
merger (jKing et al. 200^ have been advanced to inter- 



pret shGRBs. Dai et al.l (|2006[ ) argue that the final prod- 
uct of a NS-NS merger may be a heavy, differentially- 
rotating NS, whose post-merger magnetic activity would 
give rise to flares follo wing the merger eve nts. The 'off 
axis collapsar' models (iLazzati et al.ll201C)|) and the mil- 
lisecond magnetar model (Metz ger et all 12008') offer ex- 
planations for why the extended emission from short 
GRBs may resemble the prompt emission from long 
GRBs. However, it is unclear how this model can explain 
the apparent different galactic environments in which 
short and long GRBs are found. In both types of bursts, 
the photon spectrum is well fitted by the Band function, 
hinting at a kind of universality in the underlying phys- 
ical model. 
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Here, we will explore an alternative model based on 
a Quark-nova explosion occurring in a low mass X-ray 
binary (LMXB). Such an explosion can happen if the 
neutron star, in its spin-down evolution, reaches the 
quark deconfinement density and subsequently under- 
goes a phase tr ansition to the more stable strange quar k 
matter phase ()Itohl 119701: [Bodmeij [TOTlt IWittenI Il98l . 
resulting in a conversion front that propagates towards 
the surface in the detonative regime. This hypothesis 
is motivated by recent work on numerical simulations of 
the conversion of neutron m atter to strange quark mat- 
ter (jNiebergal et alJl2010bf ). Previous works that have 
focused either exclusi vely on conversion th rough strange 
matter diffusion (eg. lOlinto Ill987t iHorvat h&Bcnvenuto 
19881) or use pure hydrodynamics feg. lDrago et al.li2007l : 
Cho et ainil994[ ) do not find detonation, only deflagra- 
tion . However, in the mo st recent work on this mat- 
ter (jNiebergal et al.ll2010b[ ). we have analyzed the issue 
numerically including both reaction-diffusion and hydro- 
dynamics consistently, albeit within a ID approximation. 
We also included for the first time, the effects of neutrino 
cooling, finding numerical evidence for laminar speeds 
approaching 0.04c (much higher than previous works; c 
is the speed of light) and a possible wrinkling i nstability 
in hig her dimensions. As argued recently by iHorvathl 
(|2010i) . wrinkles cause turbulence in the conversion front 
that serve as a platform for detonation. If this is borne 
out by more sophisticated simulations, we will have found 
potentially a new engine for GRBs that is capable of ex- 
plaining several features of the GRB light curve, as we 
show in this work. This exciting possibility, combined 
with the problems faced by more conventional scenarios 
described in the previous paragraph, should make our 
model worth e xploring as a possible op tion. 

We note that lNiebergal et al.l |2010b) is already an im- 
portant step forward in resolving some of the issues with 
uncertainties in the dynamics of the Quark-Nova. The 
equation of state of quark matter, and hence the puta- 
tive quark deconfinement density is also unconstrained, 
but recentlv lOzel fc Psa ltis' (2009") have argued that ac- 
curate measurement of mass and radius for a minimum of 
three neutron stars will provide strong constraints up to 
several times nuclear saturation density. While we must 
await future observations of neutron star parameters for 
the resolution of some of these issues on the underlying 
Quark-Nova dynamics to be more confident in our model 
assumptions, we nevertheless adopt a working hypothesis 
in this paper that the detonative regime is indeed the 
likely result of a quark-hadron phase transition inside a 
neutron star. From this point on, we will follow standard 
arguments to demonstrate that a Quark-Nova occurring 
in an LMXB following the end of the first accretion phase 
leads to features (prompt emission, extended emission 
and flaring) reminiscent of those discussed above in the 
context of short GRBs as observed by Swift. The same 
engine at play inside a collapsar can naturally explain the 
many similarities between short and long duration GRBs, 
thus exhibiting universality between the two classes. In 
the end, we obtain many observed features of short GRBs 
beginning with our single hypothesis of the detonative 
transition, implying that shGRBs could well provide the 
astrophysical evidence for the Quark-nova, and that this 
mechanism deserves detailed study. 



This paper is organized as follows: In ^ we describe 
the LMXB parameters and give a brief account of the 
Quark-Nova (henceforth QN) and the relativistic ejecta 
it produces. In Sj3l we describe characteristics of the cir- 
cumbinary disk (CD). The interaction between the ejecta 
from the QN and the CD is studied in 21'^here we de- 
scribe the energetics and spectral features of the prompt 
emission. Using our model, we perform fits to observed 
GRBs that bring out the origin of the empirical Band 
function. We also derive the duration of the prompt 
emission in our model and predict an optical counter- 
part. Extended emission arising from ablating the white 
dwarf is addressed in ^ Late X-ray activity ascribed 
to the spin-down evolution of the quark star and accre- 
tion from trapped debris is explained in fj6] along with 
afterglows from the shocked CD. ^Zl contains our main 
conclusions. 

2. QUARK-NOVA IN A NS-WD LMXB 

'St aff et al.l (pOOfih . using standard equations of state 
(EOS) of neutron-rich matter, considered the likely NS 
candidates that could reach quark deconfinement den- 
sity in their core; the fiducial value was taken to be 
^ 5po, where po ~ 2.7 x 10^'*g/cc is nuclear saturation 
density. To reach ~ 5po, a NS with gravitatioiial mas s 
Mg ~ 1.8Mq (using the APR EOS: iAkmal et al.l (|1998[ )) 
is required. 

Quark deconfinement in the core of the NS can happen 
in two ways: (i) right after the neutron star is formed, if 
the above criteria is satisfied. This is not relevant to the 
model presented, but could be i mportant fo r superlumi- 
nous supernovae (|Ouved et all (I2009L \20Tdi) ): (ii) If the 
above criteria is not met, the neutron star can accrete 
mass from a companion and reach the critical mass for 
deconfinement and subsequent QN explosion. This is the 
scenario considered in our model. 

2.1. The LMXB evolution 

We b egin with the standard LMXB evolution model 
(jVerbu nt 1993), where a system with a relatively long 
orbital period (~1 day) produces a recycled pulsar when 
the subgiant progenitor of the white dwarf overflows 
its Roche lobe. In the first accretion phase, mass and 
angular momentum losses shorten the orbital period. 
Once this first accretion phase ceases, the orbital period 
continues to decay due to gravitational wave radiation. 
Eventually, the orbit decays enough for the white dwarf 
itself to overflow its Roche lobe, commencing the second, 
ultrashort orbital period (< 1 hr) LMXB phase. Here we 
are mainly concerned with the first recycling phase. 

The orbital separation of the binary, a, evolves accord- 
ing to Kepler's third law 

„ = 2.28xl0'c„,l+,)"=(^)"'(^)"\(l) 

where q = My/o/Mj^s, and Mt = Mns + AfwD the total 
mass of the system. The mass of the WD, Mwd, is con- 
strained by the expected core mass of the Roche Lobe 
filling progenitor in the LMXB phase, Mwd ^ 0.15Mq. 
We adopt Mwd ~ 0.1 Mq representing the WD's fidu- 
cial value near the end of the first accretion p hase. We 
will make use of the WD mass-radius relation (|Savoniiel 
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(fT98l ) 



^= 0.028(1 + X)V^M-/L 



(2) 



with MwD,o.i being the white dwarf mass in units of 
O.1M0. Hereafter we take X = (the Hydrogen fraction 
in the WD), since we assume a pure He WD near the 
end of the first accretion phase. 

At the end of the first recycling process (i.e. first ac- 
cretion phas e), the neutron star will reach an e quilib- 
rium period (|Bhattacharvafcvan den Heuvellll991| ) which 
is approximat ed by the Keplerian o rbital period at the 
Alfven radius ()Ghosh fc Lambl[l992l ): 



-PNS,oq. ''^ 2 msi3^g^gi?Jfg'gMMc'i^,m 



-5/7 
'NS,l-4' 



3/7 .o^ 
NS,Edd ' y"^) 



where Bns.q, ^ns,6 and AfNS,i.4 are the neutron star 
surface magnetic dipole field, radius and mass in units 
of 10^ G, 10^ cm and 1.4 solar mass respectively. The 
accretion rate onto the neutron star, ttins, is given in 
units of the Eddington limited accretion rate MEdd ^ 
10~* Mq yr~^ above which the radiation pressure gener- 
ated by accretion will stop the accretion flow. 

Neutron stars as massive as ~ 1.8Mq are not easy to 
produce in LMXBs even for initially high mass of the 
donor star, unless they we re already born as relatively 
massive neutron stars (e.g. iLin et all ()2010D ). Assuming 
up to ~ 0.2Mo can be accreted (which will spin- up the 
NS to millisecond periods)^ by the end of the first accre- 
tion phase, this puts a constraint on the minimum mass 
of the NS at birth of ~ 1.6Mq. Therefore a prolonged 
phase of accretion at the Eddington rate for > 10^ years 
is needed. Thus LMXBs with NSs born massive and 
accreting near the Eddington limit are the most likely 
systems to experience a QN explosion near the end of 
the first recycling phase. Figure [T] (upper panel) is a 
schematic representation of the LMXB system at the end 
of the first accretion phase just before the QN goes off. 

2.2. The Quark-Nova 

In the Quark- nova picture (Ouyed et al. 2002), there is 
an explosive phase transition to the more compact (u,d,s) 
quark phase, and the gravitational potential energy is 
released partly into outward propagating shock waves 
from the supersonic motion of the {u, d, s) conversion 
front. The temperature of the quark core thus formed 
rises quickly to 10-20 Me V since the co l lapse is adiabatic 
rather than isothermal (|Gentile et al.l 119931 ). As men- 
tioned in the introduction, a complete dynamical treat- 
ment of the QN, including multi-dimensional neutrino 
transport and hydrodynamics is only in the preliminary 
stages (iNieb ergal et al. 2010b) but results based on cal- 
culations incorporating the most physics suggest that a 
detonation is possible. The important outcome of the 
QN, as we now explain, is that chunks of the neutron 
star's crustal matter, rich in iron-group elements, can be 
ejected from the surface of the neutron star at relativistic 
speeds. 

^ The spin-up of a slowly rotating neutron star with 
moment of inertia /ns to the equilibrium spin pe- 
riod f'NS,cq. = 27r/f2 ~ 2 ms would require mass 
accreti on of at least ^ I-^c,Q,-^a,/ {GM-^a,R^a,)^/^ ~ 
0.2MM llBhattacharvafcvan den Heuve]||199ll ). 



2.3. The Quark-Nova ejecta 

Unlike Supernovae, neutrino-driven mass ejection in 
Quark-novae is not feasible, as neutrinos are trapped 
inside a hot and dense expand ing quark core, onc e 
it g rows to more t han ^ 2 km (jKeranen et al.l [2005f) . 
In (jKeranen et al.l I2005D . we used a typical inelas- 
tic neutrino-nucleon cross-section to show that these 
(electron-) neutrinos are absorbed in the neutron-rich 
crust of the neutron star. The corresponding diffusion 
timescale out of the conversion front is ~ 0.1 seconds. 
Although the neutrinos have a similar luminosity as in 
supernovae (using the per species number for the latter), 
there is an important difference in the outcome for mass 
ejection. Unlike the case of a supernova, where neutrinos 
can drive a wind from the surface of the proto-neutron 
star, the neutron-rich crust is much too dense for these 
neutrinos to transfer kinetic energy effectively (the grav- 
itational potential well is much deeper). Rather, they 
heat up the crust to temperatures of ~ 1 MeV. In com- 
parison to the photon-driven explosion (|Vogt et al .11 2 0041) 
which we have invoked here, energy deposition by neutri- 
nos is less than 10% of the total energy budget. There- 
fore, it is not expected to lead to much baryon-loading 
of the photon fireball. Mass ejection due to core bounce 
is also unlikely unless the quark core is very compact (1- 
2) km. A more promising alternative is mass ejection 
from an expanding thermal fireball that is a direct con- 
sequence of dense, hot quark matter pr oduced as a result 
of the conversion from neutron matter ()Vogt et all 120041 : 
lOuved et alll2010D . The fact that more than 50% of the 
gravitational and conversi on energy is releas ed as pho- 
tons is unique to the QN (jOuved et al.l 120051) . It allows 
for ejecta with kinetic energy easily exceeding 10^^ ergs. 
Depending on the conversion efficiency of photon energy 
to kinetic energy of the ejecta (the neutron star's out- 
ermost layers), up to 1O~^M0 of neutron-rich material 
can be ejected at nearly relativistic speeds. Calculations 
of mass ejection in the QN scenario accounting for en- 
ergy transfer to the c rust, give estimates o f ejected mass 
of IO-^Mq-IO-^Mq (jKeranen et al.|[2005l ). The average 
Lorentz factor of the QN ejecta is. 



rQN — 



10 



-^QN,52 



(4) 



where E^^ ;^^ 



MqnC^ Mqn-3.3 

?7o.i-£'QN,53 is the kinetic energy of the 



ejecta in units of 10^^ erg. Here 77 ^ 0.1 is the efficiency 
of energy transfer from the QN total energy (~ 10^'^ ergs 
from gravitational and phase conversion energy) to the 
ejecta's kinetic energy. The QN ejecta mass, Mqn,-3.3 
is given in units of 5 x 10"''' Mq which we adopt as our 
fiducial value. 

There are two main sources of uncertainty to calcu- 
lating the Lorentz factor (the neutrino contribution, we 
have argued above, is not more than 10% and is unlikely 
to change much even with improved calculations of neu- 
trino transport in dense matter) . The first is the thermal- 
to-kinetic convers ion efficiency rj in e q.dH), which we have 
taken as 0.1 (e.g. iFrank et al.l (|2002f )l. Fixing this num- 
ber is equivalent to fixing the amount of mass ejected. 
We are at present working on reducing this uncertainty 
by numerical modeling of the hydrodynamics of the con- 
version (from neutron to quark matter) front, including 
particle diffusion and advective forces into account (see 
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Fig. 1. — Upper Panel: A LMXB system which consists of a neutron star, a low mass white dwarf (AfwD ~ O.IMq in our model) 
and a circumbinary disk near the end of the first accretion era. Lower Panel: Near the end of the first accretion phase, the neutron 
star experiences a QN explosion ejecting the neutron star crust and leaving behind a quark star. The QN chunks are ejected during the 
QN explosion hitting the circumbinary disk material at different latitudes. Below a critical height 2cD,vap (from the midplane of the disk), 
the disk density is sufficiently high that the QN chunks are shocked above evaporation temperatures (circles with horizontal stripes). At 
heights greater than 2:cD,vap the density is low enough for the QN chunks to survive the shock and emit in the optical (circles with vertical 
stripes). A portion of the QN chunks will collide with the white dwarf causing it to partially or fully ablate, depending on the NS-WD 
separation and the white dwarf mass, when the QN explosion occurs. 

or liquid phase) The typical chunk size at birth was cal- 
culated to be Archunk,b 186 cm'0_3rb,7Pb,8 where ip is 
the iron breaking strain in units of 10^'^ (the iron break- 
ing strain is '-^ 10"'^ for solid iron and ^0.1 for the liquid 
phase), rb,7 is the radius in units of lO'' cm when break- 
ing starts and ph.s is the ejecta density in units of 10* g 
cm"^ at breakup (see §3.4 and Table 1 in lOuyed fc Leah\ 



iNiebergal et al.|[2010b[ ). Once these simulations are ex- 
tended beyond ID, a comparison of the thermal energy 
density of the conversion front versus the pressure ex- 
erted on the surrounding matter should provide a better 
constraint on rj. 

The second source of uncertainty is the EoS and com- 
position of the high density matter that is ejected. Our 
estimate for the fiducial value of ejected mass (and hence 
Fqn) is based upon a fit to the EoS of the neutron star 
crust, specifically the BP S equation of state for the outer 
crust (|Bavm et al.lll97l[) which is matched to a higher 
density equation of state for the inner crust. The equa- 
tion of state for the outer crust is reasonably robust but 
can be improved by a few percent by taking into account 
more recent models based on quantum molecular dynam- 
ical (QMD) si mulations of the nature of the crust post- 
accretion (eg. iHorowitz fc Kadaul I2009D ; astrophysical 
phenomena such as quasi-periodic oscillatio ns (QPOs) 
can also provide constraints (eg. [Watts fc Strohmaved 
|2007( ). The higher density of the inner crust is how- 
ever poorly constrained, since it involves very neutron- 
rich nuclei and possibly some "pasta" phases. We expect 
that upcoming experiments with rare-isotope beams (eg. 
FRIB; http:/ /www .jinaweb.org/ria/) will provide better 
c onstraints in t he high-density regime. 

lOuyed fc Leahy ((2009i ) studied the relativistically ex- 
panding shell of iron-rich ejecta emitted from the neutron 
star crust and found that it breaks up into numerous 
chunks because of lateral expansion forces (10"^ to 10" 
chunks depending on whether the ejecta cools in the solid 



|2009( ). The chunk's rest length l^i^^i^]^ is much larger than 
its width Archunk,b 200 cm for the solid phase (which 
we consider i n this paper) or equal to its width for the 
hquid phase (|Ouved fc Leahvl [2009| ) . In the sohd case, 
the chunks thus resemble "iron needles" moving parallel 
to their long axis. In the observer's frame the ratio of 
length to width will be contracted by I/Fqn. 

Even beyond the break-up radius, the chunks remain 
in contact with each other within the relativistically ex- 
panding ejecta as a whole. This is because the pieces 
expand in volume, filling up the space between them, 
and causing the density of each piece to continuously 
decrease, until they reach the zero pressure iron-density 
(pfc), at radius Vgcp ^ 2 x 10^ cm, at which point they 
stop expanding. From OuyedfcLeahy (2009), the typical 
chunk size at rgcp is Archunk ~ 10^ cm while its length 
is of the order of 10^ cm for the solid ejecta case. As 
the chunks continue to expand radially outwards to a ra- 
dius r > r-gop, we can associate with them a filling factor 



/< 
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We will now motivate the statistical model of ejecta 
break-up. Investigations of the fragment-size distribu- 



Quark-Novae in Low-Mass X-ray Binaries 



5 



tion (FSD) that results from dynamic brittle fragmenta- 
tion show that fragmentation is initiated by random nu- 
cleation of cracks that are uns table against side-branch 
formation (jAstrom et al.|[200^ . The initial cracks and 
side branches individually form an exponential and a 
scale invariant contribution, respectively. Both merge 
to yield the resulting fragment size distribution which is 
expressed in the general form n{s) oc s~"/i(s/si) with 
a = {2D — 1)/D {D being the Euclidean dimension of 
the space), /i a scaling function (an exponential) that 
is independent of si for f ragments sma ller than si and 
decays rapidly for s > si (jAstromI 120061 ). Here si ~ 
where A is the penetration depth of a branching crack. 
Adopting this model, the emergent picture for the rel- 
ativistically expanding shell of QN ejecta is as follows: 
cracks are initially nucleated at random and more or 
less uncorrelated locations. From these locations, cracks 
propagate in different directions. The main cracks form 
large fragments with an exponential Poisson distribution 
with a typical size determined by Archunk- The size (i.e. 
mass) of the Poisson-process fragments will be limited 
or reduced by the creation of small-size fragments (si) 
around each crack by the side-branchin g process. Th e 
final fragment size distribution becomes ()Astr6mll2006l) : 



n{s) = {l-(3) 
-l-;9exp 



s 

So, 
(si/^ 



exp 



So 



(5) 



So 



where sq — Archunk hi our model. 

Theoretical and experimental investigations 
(jSharon et all Il995t ISharon fc Fineberd [l996i ) have 
shown that a dynamic instability controls a crack's 
advance when its velocity exceeds a critical velocity 
of O.SGur where wr is the Rayleigh wave speed in the 
material (a Rayleigh wave is a surface acoustic wave that 
travels on solids). Beyond 0.36z)r, the mean acceleration 
of the crack dynamics change dramatically. At that 
point, the mean acceleration of the crack drops, the crack 
velocity sta rts to oscial late leading to the so-called YofFe 
insta bilitv (lYo fFc 1951) which leads to side-branching 
(e.g. iBuehler fc Gag ,20061) . Numerical solutions of 
the model for v > 0.36dr, exhibit the occurence of 
these branching events where the main crack sprouts 
side branches. The spacing between these branches is 
a function of the amount of dissipation. The above 
equation assumes that crack branches propagate easily 
with large penetration depth given by A = Ar 



l/D 
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Shown in Figure [2] is the FSD for two-dimensional 
fragmentation with D = 2 (i.e. a — 1.5)^. There are 
two natural scales in our model, the typical fragment 
scale Archunk (and thus mass mchunk) and, the scale 
at which the two terms/contributions in the RHS of 

^ A special case is, however, when a two-dimensional object 
(i.e. a thin plate and/or closed shells) is fragmented in a three- 
dimensional space. Models revealed that the branching-merging 
picture for two-dimensional fragmentation may be valid for the 
initial stages of fragmentation, but fragmentation may continue as 
a cascade of breakups beyond the limit of maximum fragmenta- 
tion through branching and merging of cracks. It was found that 
at later stages additiona l fragments were f ormed with a power-law 
FSD that had a ~ 1.17 l ILinna et al.|[200i) . 



Z 
Z 




0.001 



S/So 



Fig. 2. — Log-Log plot of the Quark-Nova normalized fragments 
size distribution for D = 2. The dotted top curve is for /3=0 (chunk 
distribution is purely defined by the main cracks) while the dashed 
bottom curve is for 13=1 (chunk distribution is purely defined by 
the branching cracks). The three solid curves, from top to bottom, 
are for /3 = 0.9, 0.99, 0.999 and illustrate the contribution from the 
two size distributions as described in equation JSJl. 

equation ([5]) are equal, which we call. Art. For s < Art 
FSD is dominated by the side branching. As can be 
seen in Figure [2] the transition size becomes smaller, 
i.e. TOt/rrichunk << 1 as /3 — 1. As we show later in 
this statistical model plays a key role in explaining 
the empirica l ly obs erved Band function of GRB spectra 
(iBand et all (Il99l ). 

3. THE CIRCUMBINARY DISK 

Circumbinary disks (CDs) could accompany the for- 
mation and evolution of binaries. These CDs are either 
the remains of fallback disks produce d in the supernovae 
that formed the compact object (e.g. IWang et al. Il2006t 
ICordes fc ShannonI l2008l ) . or material injected into cir- 
cumbinary orbits during the pr ocess of mass loss by the 
Roche lobe-filling companions (|Dubus et al.ll2002[ ). The 
possibility that the CB disk is the remnant of the com- 
mon envelope evolution phase of binaries cannot be ruled 
out either. Due to the tiny CD mass that would re- 
sult, even after billions of years, such disks would be 
difficult to detect in the optical and infrared. Nev- 
ertheless, the detection of excess mid-infrared flux in 
few LMXBs in quiescence suggests that dusty circumbi- 
nary material might be present around some L MXBs 
(|Muno fc Mauerhanl2006[ l (see also lBOTdeJ (|2010[ )'). This 
remains to be confirmed. Here we consider (and assume) 
the case in which the binary system is surrounded by 
stable circumbinary disk/material continuously replen- 
ished during the process of (near-Eddington) accretion. 
As we have said, building massive NSs in LMXBs would 
require near Eddington accretion rates. This, we argue, 
could favor the formation and the sustainment of the CD. 

There are different possible configurations for the CD 
such as a flat or flared disk. Here, we assume the re- 
plenished disk to be a geometrically thin, viscous Kep- 
lerian circumbinary disk. Is scal e height, H\r,, at t he in- 
ner edge, Rin, is Fin ~ 0.03i?in (|Belle et al.l 120041 ). The 
disks could lie as close to the center of mass of the bina- 
ries as i?in 1.7a, at which point they wou ld be tidally 
truncated ()Taam et al.l 120031: iDubus et al.l I2004D . The 
dependence of the vertically integrated surface mass den- 
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sity, E on radius R, is taken from (|Taam fc Spruitl[2001[ ). 
particularly Figure 8 in that paper, from which one can 
approximate E ~ Ein{Rin/R)- The mass of the disk is 
then 

/•flout 

McD = / ^{R)RdR ^ Ei,i?2ji?„^j/i?i„ _ 1) (6) 
"'-Riu 

-1.5 X lO"^M0Ei„,5a^oCcD,out 

where CcD,out = (-Rout,i2/^^in40 - 1) • For a disk that is 
vertically isothermal, at a given radius i?, the equation of 
hydrostatic e quilibrium combi ned with the radial profile 
above yields (|Belle et al.ll2004D : 



PCD = 



Rjn 

~~R 



cxp 



(7) 



Unlike a Shakura-Sunyaev disk, in CDs, the opening 
angle, H/R, is expected to be roughly constant (here 
— 0.03). Recalling that i?in ~ 1.7a, the above then yields 



PCD ^ Pin 



exp 



(8) 



where pin 10 g cm ^ (Ein^s/aio) and Ein^s = 
Ein/10^ g cm^. The column density as seen at a viewing 
angle 90° — i = ^view = Zview/R, where i is the inclination 
angle of the system is 



1 



' 6 X lO^^Es exp 



PcvdR 



(9) 



n2 
view 

26*? 



where niH is the proton mass and /icD ^ 1.2 is the 
mean molecular weight for the CD material which we 
assume has a solar composition. 

4. PROMPT EMISSION 

We first address the fate of the QN chunks when 
they impact the CD. The chunks shock the CD mate- 
rial, and in turn, undergo a reverse shock that heats 
them to a temperature found from pressure balance 
-Pchunk = ^CD,sh.- The prcssurc in the shocked C D ma- 
terial, PcD.sh.i is fou nd from the jump conditions (|Russol 
ll988LIIwamotolll989l) . We get, at i? = i?in, 



fceTchunk ^ 30 keV 



V r2 

^in,5J- QN,10 



exp 



2m 



(10) 



Noting that non-degenerate ir on will vaporize if heated 
to > 0.3 eV (jCRC-PressI 120051 for vaporization temper- 
ature of iron at normal density), the rapid fall-off of 
the CD density given by eq.® implies that the vapor- 
ization is already efficently accomplished by interaction 
with the inner edge of the disk. At the inner edge of 
the CD, we can define a critical disk height, 2cD,vap. 
(and corresponding density pcD,vap.), below which the 
chunks will start vaporizing upon impact (see Fig. [T] for 
a schematic illustration of the process). This critical 
height is zcD.vap. ~ 4.8-ffin with a corresponding CD 
density pcD,vap. ^ 5.8 x 10^^ g cm~^. It follows that 



the solid angle extended by the critical surface defined 
by p > PcD,vap. is then 



27ri?in22:cD. 



vap. 



'iirRi 



Ru 



0.15 



(11) 



The amount of QN ejecta material contained within flc 
is Mqn.c = f^cMQN ^ 7.5 X 1O-5M0Mqn,-3.3- 

The time it takes a given chunk to be completely 
eroded is 



1/2 



chunk 



4.3 X 10-^ s 



QN 



Vl/2p3 
^^in^S^ QN40 



(12) 



where Vs = a/ (fcBT'c/2mH) is the speed at which the heat 
wave crosses the chunk. The chunk's Thompson opti- 
cal depth, ^ l/(?T-chunkO'T), is much less than the typical 
chunk size, Archunk 10^ cm (the chunk's density is 
about ~ 10^ g cm~^ by the time they reach the CD; see 
§5). During their extremely brief life as solid entities in 
the CD, most of the chunks will emit as blackbodies at a 

rate Lchunk = ^chunkO-T^unk where Achunk ^ '^^r'^hnni,- 

For a typical chunk, ichunk ~ 10^° erg s~^ T'c',30 where 
the chunk's temperature is given in units of 30 keV. The 
total luminosity in blackbody emission can be estimated 
by recalling that the total area extended by the chunks is 
- Ana^ so that Lbb nATra^crT* - lO*^ erg s-^r^gg. 
The corresponding isotropic blackbody contribution from 
all the chunks is Lbbjso. — Lbb/^c ~ 10^" erg s~^. An 
observer would see a blackbody peaking at FqNTchunk- 

A given chunk will continue to plow CD material even 
after complete vaporization and erosion since the va- 
porized material will retain the directional motion and 
size distribution of the incoming chunks. The vaporized 
chunk material will eventually enter a decelerati on phase 
once it plows a mass of ^ m-chunk/FqN (e.g. iRhoadsl 
Il997f ) in CD material. This puts a constraint on the 
minimum mass of the CD disk necessary to stop the 
QN ejecta within flc- It is Mco.min = -Mqn.c/Fqn ^ 
7.5 X 1O~^M0Mqn,-3.3. In our model, this translates 
to a minimum disk size i?out,inin = ^stop given by 

CcD,stop — -Rstop/^in ~ 500 Mqn, -3. 3/(^1,1, 5aiorQN,lo)- 

This critical radius is important for the kinematics and 
dynamics of a given chunk in the CD. Assuming a CD 
with i?out > ^stop implies that most of the vaporized QN 
ejecta will lose its momentum and kinetic energy to the 
CD material before it reaches the outer edge of the disk 
(see S21). 

4.1. Spectrum: Origin of the Band Function 

For the prompt emission, we isolate two broad regimes 
(in addition to the blackbody emission from the heated 
chunks) in our model depending on whether the shock 
temperature of the chunks exceeds the nuclear dissocia- 
tion temperature (binding energy) for iron or not. 

4.1.1. Regime 1 - MeV photons, transient Fe lines 

Upon impact, the chunks are heated to temperatures 
above 0.3 eV but below 8.8 MeV, in which case they 
vaporize without dissociation. Prompt emission in our 
model is tied to the shocked CD material, whose tem- 
perature TcD.s can be estimated as follows: The chunk 
plows through the CD while vaporizing into a hot stream 
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of gaseous particles which retain the directional motion 
of the incoming chunks. In this case, in the observer's 
frame, 



PCD ^ 2 dMcD 

—7- i QNWchunkC ~ 

Pchunk AiCDTOH 



-fcBlcD.s • 



(13) 



We have dMcD,piow=^chunkPcD'i-RcD,obs. where 
^chunk is the area of the chunk and dRco.ohs. is the 
distance traveled by the vaporized chunk as measured 
by an observer. The above equation is valid as long as 
dMcD.piow < "T-chunk/rQN- The chunk's rest mass is 

Wchunk = Prhiitik^cliunk^'hiink with p^hunj^ — TgNPchunk- 

Replacing in egn. dlSt yields 



knT, 



bJcd,s ^ fJ-CDrriHC 



chu 



CD,obs 



A limit on the observed temperature is then. 



fceTcD.s > lOkeV^icD x 



chunk, 5 
'^-RcD,obs,10 



(14) 



(15) 



where c?i?cD,obs ~ c x ~ 10^° cm to account for a 
typical GRB duration dicRB ~ 0.3 s (see ^4.3p . The 
temperature given above is the peak temperature after 
the chunk has plowed through most of the CD. Higher 
peak temperatures occur in the early stages of the plow- 
ing process. 

Ean. (|14l) is illuminating - it implies a direct propor- 
tion between the temperature of the shocked material 
and the linear size of the chunk. Since the linear size is 
distributed according to Eqn. ([5]) , the temperature of the 
shocked CD material will also display the same distribu- 
tion. 

It is interesting to note that several GRBs show spec- 
tra (more precisely, vF(i /) curves) that ap pear to follow a 
log-normal distribution (|Band et al.lll993[ ). According to 
iBrownfcWohletd (|1995| ). fragmentation of a massive ob- 
ject that proceeds by formation and branching of cracks 
results in a distribution of masses (per unit lo garithm 
in m ass) that follows the WeibuU distribution (jWeibullI 
119391) , which looks almost identical to the log- normal dis- 
tributi on. The WeibuU distrib ution in mass can be writ- 
ten as (jBrownfe Wohletzl 1 1 9951 ) 



5+2 

m''n(m) — Nmi I — ) exp 
\ mi ' 



(m/mi) 



1 



(16) 



where S — —D/3 is a fractal dimension and mi a mass 
related to the average fragment mass in the distribution. 
The peak of the distribution (called the most proba- 
ble mass; in our case mchunk) is given by mpcak/"ii — 

(2 -I- (5)i+^; TV is a normalization factor related to the 
total number of fragments. Note the formal similarity 
to the first term in eqn.®, derived by Astrom (2006), 
who based his expression on a more complicated frag- 
mentation situation involving the competition between 
side-branching and cracking, which leads to the second 
term in eqn.®. 

Therefore, we put forward the intriguing suggestion 
that the GRB spectrum arises due to photons from the 
shock-heated CD material whose temperature distribu- 
tion is inherited from a fragmentation distribution of the 



QN ejecta. We assume a one-to-one relationship be- 
tween the temperature of the shocked CD material and 
the resulting radiation energy; in this case, the prompt 
photon energy Ej will inherit the fragmentation distri- 
bution. Although a full calculation of the non-thermal 
spectrum of the photons is beyond the scope of this 
work, we test this idea by performing fits to some ob- 
served GRB spectra by assuming that the photon lumi- 
nosity {vF{v) = E^n{E^)) follows the WeibuU distribu- 
tion given by (fT6|) . 



E-y 



(5+2 



exp 



<5+n 



S + 1 



(17) 



where is a normalization factor and Eq the average 



photon energy in the distribution, 
tribution occurs at 



The peak of the dis- 



^pcak= (2 + <5)T+7 £;o . (18) 

In our model, the assumed one-to-one correspondence 
between Tcd.s and E-^ implies that -Epcak evolution in 



time is given by eq.dlll). The fits to GRB 000218, GRB 
000727, and GRB 980706A are shown in Figures 
In reality the spectrum evolution will be related not 
just to the chunk distribution but also the vaporitza- 
tion rate. We expect the softening of the spectrum to be 
caused mostly by the erosion of the fast moving chunks 
which will be heated the most and thus radiate the most. 
The total number of chunks should also decrease in time 
with smaller chunks most likely disappearing (i.e. fully 
eroded) first. The case of GRB 000727 is an interest- 
ing one since its observed Eq increases in time. However 
note the data is for shorter integration times than for the 
other two candidates (see Table 1). In our model, there 
is a minimum accumulation time (presumably of the or- 
der of milliseconds) before the full mass distribution (i.e. 
the Band spectrum) builds up. The excellent fits suggest 
that the underlying physical picture relating the spec- 
trum to the ejecta fragment size or mass distribution is 
probably correct. In reality, the fragmentation process 
may be more complicated (e.g. eq.®), but the WeibuU 
distribution appears to be an excellent approximation. 

There are specific predictions for this component 
within our model: (i) The peak photon energy (related 
to ^JijujjIj.) should evolve in time according to equation 
(|14p ; (ii) The prompt emission is composed of contribu- 
tions from ^ 10^-10^^ "blobs" of emitting CD material 
plowed and heated by the relativistic vaporized chunks 
material. A given pulse in our model consists of different 
"blobls" of heated CD material emitting simultaneously. 
Thus a pulse should also carry the signature of the dis- 
tribution in mass of the chunks (i.e. the Band function); 
(iii) finally, the prompt emission (regime 1) should have 
signatures of iron-group elements (Fe, Co, Ni) brought 
in by the chunks, but this should be transient since the 
vaporized material will cool rapidly. 

4.1.2. Regime 2 - GeV photons, no Fe lines 

Chunks heated above 8.8 MeV will be dissociated into 
nucleons and lose their identity by mixing with the CD 
matter. From Eqn. ()10p . this regime could occur for sys- 
tems with Fqn > 100 (i.e. with A/qn < 5 x 10~'^Mq). 

The mixed CD and QN material will have a common 
temperature and will have some neutron excess (since 
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Fig. 3. — Model fit (soUd curves) to observed photon spectrum {vF^) 
with data from Mazcts ct ah (2004) for GRB 000218 (top panel), GRB 
000727 (middle panel), and GRB 980706A (bottom panel). The top 
and bottom data points, in each panel, are for different integration 
time (see Table 1). 

the QN ejecta inherits its composition from the neutron 
star crust) with a mean molecular weight /i ^ 1. We 
expect a given chunk to instantly transfer a significant 
fraction of its kinetic energy into heating m.chunk/rQN of 
CD material. The resulting thermal energy per nucleon 
is then 

fcBTnuci. - 100 GeV r^j,^o . (19) 



For such energies which exceed the threshold for pion 



Fig. 4. — Model fit (solid curves) to observed photon flux n{E-y) 
with data from Mazcts et al. (2004) for GRB 000218 (top panel), GRB 
000727 (middle panel), and GRB 980706A (bottom panel). The top 
and bottom data points, in each panel, are for different integration 
time (see Table 1). 

production, processes such asp + p— )-p+p + 7r°,ri+p— > 
p + p + T:~,p + p—^p + n + 7r+ can occur. Charged pion 
decay produces electrons and positrons which can anni- 
hilate to form photons, and neutron pions can produce 
photons by 7r° — 27. The luminosity can be estimated 
from 

(E-^) ^'^^""^^'^P 27ri^iniJcD^chunkCCrpp (20) 
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TABLE 1 

Model fits to observed photon spectrum (GRB 000218, 
GRB 000727, and GRB 980706A) using eq.|[T3. The fits are 

PERFORMED FOR OBSERVATIONS AT TWO DIFFERENT INTEGRATION 
time, tint.- 



GRB 


tint. 


<5 


Eo (keV) 


-Epcak/ Eo 




000218 


0-256 ms 


-0.2 


520 


2.1 


0.164 




256 ms-5.888 s 


-0.5 


420 


2.3 


0.027 


000727 


0-128 ms 


-0.33 


73 


2.2 


3.2 




256 ms-768 ms 


-0.65 


80 


2.4 


5.0 


980706A 


0-256 ms 


-0.39 


760 


2.2 


0.2 




256 ms-8.448 s 


-0.67 


580 


2.4 


0.009 



--2 X 10 ergs/s I]i„,5/chunk,5aio 

where we have assumed a typical cross-section lor pp col- 
lisions in the te ns ol GeV range of <7 p„ ~ 10~^^ cm^ (i.e. 
100 mihibarns, iBlattnig et alTl2000D and a typical pion 
energy (i?^) ~ 150 MeV (threshold production). 

We offer specific predictions for this high-energy com- 
ponent within our model: (i) We expect this component 
also to show traces of the chunk size distribution and 
partly follow a log-normal-like distribution. However, 
chunks heated to extreme temperatures (>> 1 GeV) 
will start to move away from the band function; (ii) it 
should have no signatures of Iron-group elements since 
the ejecta is dissociated completely to nucleons; (iii) GeV 
and higher energy photons should be observed only in 
systems with Fqn > 100. 

For the two regimes isolated above, the CD mate- 
rial between photons and the observer must be optically 
thin. The radiation will scatter most off of free elec- 
trons in the ionized CD material. The Compton opti- 
cal depth is Ncbo-c ~ 2.7 x E5 exp (-01-^^/20'^) where 
<Tc — (3/8)(7t In (a;)/x^ is the Compton cross-section, 
(7t — 6.652 X 10~^^ cm^ the Thompson cross-section, 
and X is the ratio between the prompt photon energy and 
the electron rest-mass which is ~ 1 GeV/(mcC^) >> 1. 
The shocked matter is therefore Compton-optically thin 
to the high-energy photons. For the low energy pho- 
ton regime, the Thomson optical depth is too high and 
this should negatively affect the ability to see the shocked 
chunks or the shocked disk gas since the deposited energy 
should thcrmalizc over the whole disk. The disk geom- 
etry and vertical scale height, in our model, is assumed 
to be the simplest with H/R = constant. However other 
disk configurations are possible, such CD with a 

flared surface; H cx R^^^. This is a great advantage 
since a large solid angle is available for the chunks to be 
able to plow at an angle into increasing smaller density 
and thus be able to radiate in the optically thin regime. 
Specifically, for such a flared disk/atmosphere most of the 
chunks would be flying through the larger volume of the 
disk/atmosphere so most energy would be deposited in 
the region that is dense enough to significantly slow the 
chunks. Up high in the flared disk/atmosphere region the 
density is so low that the chunks fly through unimpeded. 
In summary, a -ff cx i?^/^ would make a large visible im- 
pact surface available with a a longer lived black body 
from the optically thick part and a short lived brighter 
part from flared region of the disk. 

4.2. Energetics 



The energy deposited by the QN ejecta (into ~ 
^^QN,c/Fqn of CD material) gives us an estimate of the 
prompt GRB emission. The total thermal energy pro- 
duced which we assume to be roughly equal to the to- 
tal GRB energy is E'grb ^ Fqn x (Mqn,c/Fqn)c^ 
1.5 X 10^" erg Mqn,-3.3 where we have used AIqnm = 
rJc-MQN- The effective isotropic energy in our model is 
thus £'GRB,iso. ^ Eghb/^c ~ 10^^ erg Mqn,-3.3- Ob- 
served short GRBs also have i?GRB iso up to lO'^^erg 
(lNakaHl2007[ ). 

In our model, the disk only subtends an angle ^ 15% 
of the sky. Thus, at most ~ 10^^ ergs of the QN kinetic 
energy 10^^ erg intersects the disk to be radiated. It was 
shown by (Zhang et a l. 2001) that the radiative efficiency 
{Ej/ {Ej + Ek)) could vary from a few percents to more 
than 90% depending on how the kinetic energy is calcu- 
lated. These studies are based on key assumptions about 
(i) the unknown shock parameters; (ii) the required de- 
tailed afterglow modeling with its own limitation and 
(iii) assumes the fraction of the electron energy in the 
internal energy of the shock to be 10%. Nevertheless let 
us assume an efficiency of ^ 0.1 for short GRBs which 
is also the number derived in an independent study by 
(Berger 2007). A 10% radiative efficiency would imply 
an E^ ^ lO'"^" erg in true energy release corresponding 
to an isotropic energy release of E^^iso. ~ 10^^ erg. This 
number could be higher for disks with higher solid an- 
gle (e.g. flared CDs) and/or for QN ejecta with higher 
kinetic energy. 

4.3. Duration 

In our model, emission duration is related to the time 
it takes the QN ejecta to be stopped by its interaction 
with the CD: 

_ (i?stop_-_Rin) 1 ^ „ MqN-3.3 . 

IGRB ^ ^-^2 ^ S t;^ • (21) 

^Ci QN ^in,.5aioi QN,io 

The time-variability can vary from microseconds to a 
fraction of a second depending on the number of chunks 
(out of the ~ 10^-10^^) that emit simultaneously at a 
given time. 

4.4. The optical counterpart 

The QN chunks hitting the CD in regions where pcD < 
PCD.vap. (i-e. at heights z > 2;cD,vap.; see Fig. [T]) wiU 
not be fully vaporized, nor do they expand significantly. 
Rather, they pass through the CD, effectively punctur- 
ing it. During this interaction the temperature of a QN 
chunk is determined by shock heating, and will not ex- 
ceed ~ 3 eV Fqn, 10 as seen by an observer; thus emis- 
sion would be in the optical band. The observed opti- 
cal emission is related to shock efficiency with emitted 
energy Eopt. ^ (0.3 eV x MQN,opt./(A'cm-H)) x F^^ ^ 
10^^ ergs Fqp^ ]^QAfQN,-3.3 where /lo = 2 is the mean 
weight per electron and MQN,opt. = f^Opt.-^QN with 
f^Opt. = (^;cD,Opt. - ^:cD,vap.)/^in ^ 0.3i7i„/i?in with the 
maximum height zcD,Opt. is calculated from ea. ([TU|) us- 
ing the lower limit of the optical band. Our model pre- 
dicts optical emission to appear simultaneously with the 
prompt emission since the chunks penetrate all regions of 
the disk at the same time. The optical and prompt emis- 
sion are both composed of short pulses from the impact 
of the chunks striking the CD. 
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4.5. Afterglow activity 

The innermost CD material acts as a buffer for the 
dissociated and pulverized chunks. The cross-sectional 
area of the chunks increases rapidly sweeping up more 
CD mass, which leads to a run-away process. As they 
expand, they start increasing the covering factor of the 
blast. If enough chunks get dissociated a covering/filling 
factor of unity would eventually be reached. Energy- 
Momentum conservation implies that the mixed ejecta is 
ejected with speed, V[, given by 



1 



r2 - 1 

^ QN ^ 

I AfcD 



M, 



QN,c 



M, 



QN,c 



-r, 



CD 



QN, 



(22) 



where f3[ — Vf/c. Depending on Mqn.cA'^cd, the 
mixed ejecta resembles a massive shell which expands 
at mildly relativistic speeds. External shocks with the 
circumstellar matter will slow down further the mixed 
ejecta and should produce an afterglow, which could 
lasts for months. This is similar to afterglow activity in 
the internal-external shock model of GRBs (Piran 2001) 
except that the prompt emission is caused by the chunks 
rather than by internal shocks. Additional afterglow 
activity can also occur if the non-dissociated chunks 
interact with matter beyond the outer edge of the CD. 



To summarize, the QN naturally provides the rel- 
ativistic "bullets " req uired by so me GRB models 
(|Heinz&:Begelmanlll999l:IUmedall2000[ ). However, the sur- 
rounding environment and emission region is very differ- 
ent in our model. The interaction between the relativistic 
QN chunks and the CD (as well as the chunks' composi- 
tion) is the key to the successes of our model - the origin 
of the Band function, as well as energetics, variability and 
duration which are consistent with observations. The in- 
teraction between the QN chunks and the WD is another 
feature of our model which we explore next. 

5. EXTENDED EMISSION 

Thus far, we have focused on the interaction between 
the chunks of ejecta and the CD as a means to explain 
the prompt emission. We now investigate the interac- 
tion between the chunks and the WD. As we show below, 
the WD can be ablated by the highly energetic chunks, 
leading to properties reminiscent of the extended emis- 
sion (EE) observed in some shCRBs. The density of the 
QN ejecta at a distance r from the point of explosion is 
given by eq.(B3) in Ouyed&Leahy (2009). A combina- 
tion of mass conservation and thermal spreading of the 
QN ejecta thickness, Ar, gives pqn ~ 10'^ g cm~^ po,i4 ^ 
Aro.4/ [t^q'^Mq^ 3 3) where the density at ejection point 

(the surface of the NS), po, is given in uni ts of 10^'^ g 
cm~^ (for an Mqn '-^ 5 x 10"'* ejecta; iDatta et al.l 
II995I ). The ejecta thickness at ejection is Atq ~ 10"* cm 
or Aro,4 ~ 1 in units of 10** cm. The chunks' density at 
the WD location (r = a) is then ~ 10'^ g cm~^ po,i4 x 
Aro,4/(a^Q'^MQ^ _^ 3) which is of the order of the WD 
density pwD ~ 6.3 x 10^ g cm~^ ^^wd i- We thus 
expect the QN induced shock to pass through the 
WD rather than flowing around since the "crushing 



time" ~ (Rwu/c) x (pwd/pqn) ~ {Rwu/c) << 1 s 
(jKlein et al.lll994[ ) . Thus it is reasonable to assume that 
all of the intercepted kinetic energy will go into heating, 
rather than shocking on and flowing around, the WD. 



5.1. 



Energetics 

The total thermal energy deposited by the chunks im- 
pacting the WD is fiwD-Eqi^, given by 



£;th ~ 9.3 X 10^9 erg E^^^,,M^'^%,aY,^ 



(23) 



where f^wD ~ RwdH^^^) ^he sohd angle subtended 
by the WD. If ii^th exceeds the gravitational binding en- 
ergy of the WD 



1.4 X 104« erg M^^^ (24) 



then the WD would be ablated as a result of the QN 
explosion. The condition above imposes a nearness re- 
striction on the WD as 



a < oth. 



8.2 X 10*° 



(-^01^52) 



1/2 



3/2 

WD, 0.1 



(25) 



For systems with a < ath., one should observe both 
prompt and extended emission. 

The temperature per nucleon of the ablated WD will 
be an important quantity that determines the spectrum 
of the extended emission. It can be estimated as 



B^WD 



133 keV 



MWD,4/3-t^QN,52 

2 n f^/^ 
^10-'"WD,0.1 



(26) 



where /xwd is the mean molecular weight in units of 4/3 
(for a pure Helium WD). This implies that WD material 
is ejected at speeds 



Vt 



th. ,ejec. 



3 X 10-^ km s" 



-'^QN,52 



1/2 



aio^^wao.i 



(27) 



5.2. Spectrum 

Nuclear burning timescales are much shorter than dy- 
namical timescales. Depending on Twd, various nuclear 
burning processes are expected to take place. As for the 
case of prompt emission, we can isolate the following in- 
teresting regimes: 

5.2.1. Regime 1: T^^d > 0.26 MeV, Ni production 

From equation (j26p . this temperature regime corre- 
sponds to 



aio < am = 0.7 



,1/2 pKE 1/2 

^^WD,4/3^QN,52 



5/6 

WD, 0.1 



(28) 



In this range, the WD temperature Tvvd > ?Si ^ 3 x 
10^ K, the temperature for Silicon burning. As a result, 
a-burning can occur all the way up to Nickel-56. The 
sequence of a-burning from Helium through to Nickel will 
release on average 7-8 MeV per nucleon, resulting in an 
energy release ~ 2 x 10^° x Mwd,o.i erg, which is enough 
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for complete disruption of the WD. The processed nuclei 
are ejected at speeds Viiuc.,cjcc. ^ 

2.1 X 10"^ km s"\ The 
high speeds and nuclear processing imply that one should 
observe in this regime 

• Broad emission lines ('-^ Kiuc.,ojcc.) early in the 
spectrum. 

• No Hydrogen or Helium lines, weak or no Silicon 
line. 

• Calcium, Nickel and Iron group lines 

The signature of a QN in such a binary system is then 
similar in some respects to what is observed in Type Ic 
supernovae. 

5.2.2. Regime 2: 50 kcV < Two < 0.26 MeV, Silicon 
production 

From equation (|26|) . this condition corresponds to 



am < flio < asi = 1.6 



/^WD.4-'^i 



QN,52 



1/2 



5/6 

WD, 0.1 



(29) 



This is sufhcient for Carbon-burning and Oxygen burn- 
ing, but not enough for Silicon burning. The observa- 
tional signature of a QN in such a system would be 

• No Hydrogen or Helium lines 

• Oxygen, Magnesium lines 

• Strong Silicon line if Oxygen burning also occurs 
(T > IQ'^K) 

This is similar in some respects to Type la SNe except 
that no Nickel can be produced. However, the QN ejecta, 
if not completely destroyed on impact, will still display 
weak Fe group lines since it is Iron-rich. Thus, these lines 
can show up in the late-time spectra making it closely 
resemble the characteristics of a type la SN. 

5.2.3. Regime 3: 8.6 keV < Twd < 50 keV, Oxygen 
production 

This is sufhcient for Helium ignition but not for Carbon 
burning. From equation (1261) . this regime corresponds to 



asi < aio < ao = 3.9 



f^WD,4-^QN,52 



5/6 

WD, 0.1 



(30) 



The total nuclear energy released ~ 1.8 x lO^*' x 
-^WD.o.i erg is once again enough for complete disrup- 
tion of the WD, but the processed nuclei are ejected at 
slower speeds Kmcojcc. ^ 1-3 x 10** km s~^. a-burning is 
the most important process here, so that nuclei 0,Ne,Mg 
are likely to be synthesized. Since Silicon burning cannot 
happen, elements such as Calcium cannot be synthesized 
in such an explosion, but can be brought in by the QN 
ejecta since r-processing in the decompressing ejecta has 
been shown to prod uce abundance peaks around Cal- 
cium and Titanium (jJaikumar et al.ll2007[ ). We expect 
signatures of such a binary to be 



• No Hydrogen, no Silicon but weak Helium lines 
(leftover Helium). 

• Oxygen, Magnesium lines and perhaps Calcium 
lines 

This is similar to the spectrum of some Type Ic super - 
novae that show weak He lines (Mathcso n et al.ll2000l) . 
An important point in our model is that Calcium and 
Oxygen originate from different sources: Calcium from 
slow-moving QN ejecta that do not get completely pul- 
verized upon impact and Oxygen from Helium-burning 
of the white dwarf. 

5.2.4. Regime 4-' Two < 8.6 keV, thermal ablation 

Nuclear burning is not possible in this system since He- 
lium ignition cannot happen in this temperature regime, 
which corresponds to 

aio > ao (31) 

If a > flth, the extended emission will be weak, but if 
a < flth, thermal ablation can still take place, according 
to the energetic arguments discussed in 15.11 We expect 
signatures of a QN in such a binary to be 

• Strong Helium lines, no Hydrogen or Silicon lines 

• ck-elements such as Calcium and Titanium 

This is similar to the spectrum of Type lb supernovae, 
except that no Oxygen line should be seen. We note 
that the distinctions between the regimes may not be 
as clear cut as presented here, rather there should be 
a continuum of possibilities. A detailed calculation of 
the observed spectral features is beyond the scope of this 
work, but will be the subject of further studies. 

5.3. Duration and variability 

In our model, the extended emission immediately fol- 
lowing the prompt emission (in the allowed regime) 
consists of emission of softer X-ray photons (Two < 
?CD,obs.) with total energy output ~ 10^" erg. Whether 
the WD is ablated thermally or ignited via nuclear re- 
actions, the energy will be released within a dynamical 

time scale Td ~ (Gpwd)~^^^ ~ 46.5 s/Mwd,o.i- 

The extended emission associated with short GRB tails 
are highly variable and in some cases on timescales less 
than ^ 1 s much less than typ ical duration (~ 100 s) of 
the tail feg. lNorris fc Bonneljl2006( ) . Currently, it is not 
clear how such variability can be accounted for in our 
model. One possible future avenue is to explore shock 
driven instabilities during the "crushing" of the WD or 
mechanisms that could drive radial oscillations of the 
WD between its RL radius, -Rwd.rl, and i?wD- 

6. LATE X-RAY ACTIVITY 

6.1. Nickel decay 

Systems experiencing a QN explosion while a < ONi 

will process most of the WD into Nickel, which will decay 

^^Ni — Co through ener getic photons with energy ~ 

0.8 MeV (jPinto et al.ll2001l ). The resulting energy release 

is ^ 1.7 X 10"*^ X MwD.o.i erg in the X-ray. If released in 

10^ s (rNi ~ 6 days) this corresponds to a luminosity of 
- 10^3 gj.g g-i. 
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6.2. Accretion onto the Quark star 

The fate of the WD ablated material depends on 
whether it can escape the system. Comparing the escape 

speed, Vcsc. 2 x 10^ km s"^ -^Qs^i.s/oio^' to Vth.,ojec. 
and Vnuc.cjcc. one can argue that, in general, part of the 
WD ablated material can be trapped. The trapped ma- 
terial could then be accreted onto the QS which could 
lead to X-ray activity beyond the EE discussed above. 
The corresponding duration can roughly be estimated as 

the free-fall time, ~ 50 s x a^ii^ /Mq^^ g, where the 
QS mass is given in units of 1.5Mq. If enough material 
is accreted in the process, the QS might turn into a black 
hole. The conversion of the QS to a black hole should 
lead to bursting events (X-ray spikes) occurring shortly 
after the EE. For systems with a < cni (with signatures 
of Nickel burning), the free-fall timescale and therefore 
the duration of the accretion activity from accretion is 
brief (< 50 s). They should also show a late-time (~ 10^ 
s) Nickel decay bump lasting for tni ^ 6 days. On the 
other hand, for a > ONi, the duration of the late X-ray 
activity could last as long as ~ 10^ seconds and no ^^Co 
or ^^Fe lines should be present. 

6.3. Quark-star spin-down 

In the event that the QS does not turn into a black hole, 
it becomes a source of late X-ray activity. The QN com- 
pact remnant is a QS which consists of a vortex where 
the interior magnetic field is confined (jOuved et al.l 
|2004|) . These vortices are expelled as the star spins 
down releasin g and dissipating its int erior magnetic field 
(|Ouved et al.l l2004t iNiebergal et al.l l2010a^ ■ For opti- 
mum conversion of spin-down energy to radiation, the 
resulting X-ray luminosity is 



Lsd--3.75 X lO^s erg s'^ 

2 



(32) 



Bo 



IQis G 



2 ms 



1 + - 

T 



-5/3 



where Pq = Pns.oq. (see eq.Q) and Bq are the birth 
spin period and magnetic field strength of the QS respec- 
tively; the QS magnetic field is given in units of 10^^ G as 
found from studies of mag netic field generation in quark 
matter fe.g. llwazakill2005l). The characteri stic spin-down 
time (in seconds) is ( Niebergal et al.ll2006[ ). 



r., = 3.3x103s fl^V m (^] fl^y 
\ Bo J \2msJ \1.5MqJ \Rqs ) 

(33) 

The curve from spin-down is flat for up to Tsd before 
decaying at a rate t~^^^ (different from the usual t~^ for 
the neutron star case because vortex expulsion in the QS 
case leads to faster magnetic field decay). The flat seg- 
ment could end abruptly in the event the QS collapses 
to a BH because of spin- down and the su bsequent in- 
crease of its core density (jStaff et al.ll2006D . Finally, we 
note that the spin down energy will be released prefer- 
ably along the equator and should naturally be observed 
together with the prompt emission. 



In this paper we presented the idea of a Quark- Nova oc- 
curring in an LMXB (NS-WD binary) during or near the 
end of the first accretion phase. Our basic assumption 
of a detonative phase transition occurring in the massive 
neutron star in an LMXB leads to many interesting con- 
sequences. We found that several features of the short- 
hard GRBs can arise naturally in our model, with the bi- 
nary separation playing an important role. Prompt emis- 
sion occurs due to the interaction between the Quark- 
Nova ejecta and the circumbinary disk. Adopting a size 
and speed distribution for the chunks can explain the 
range of energetic photons observed (keV-TeV) as well 
as the duration of the prompt component. Most im- 
portantly, a statistical approach to the ejecta break-up 
and collision with the circumbinary disk could in princi- 
ple lead to a light curve closely resembling the empirical 
Band function. 

Extended emission in our model comes from ablation 
of the WD, either thermally or via nuclear ignition. In- 
terestingly, the spectral features depend principally on 
the binary separation, and are expected to display many 
similarities to type I SNe. Thus, our model offers a plau- 
sible unification all the Type I Supernovae. Late X-ray 
activity is explained in terms of Nickel decay and accre- 
tion of gravitationally trapped material onto the Quark 
star. It may or may not become a black hole, and X- 
ray activity can result in either case; in the former it 
will consist of X-ray spikes while in the latter, it will be 
the quiescent spin-down luminosity of the QS. Afterglows 
are explained in terms of synchrotron radiation from the 
rapidly expanding shocked circumbinary disk material. 
In all cases, the energetics provide the correct order of 
magnitude. 

The rarity of short GRBs calls for rare situa- 
tions/events. Our model suggests Quark-Novae in 
LMXBs born with heavy Neutron Stars are the likely 
engines of the burst. The universality of our model can 
be confirmed if it turns out that indeed heavy neutron 
stars (> I.GMq) in LMXBs are associated with or favor 
the formation of circumbinary disks. The characteristics 
of the disk (e.g. size, density etc ...) are expected to vary 
slightly from one system to another which could explain 
differences in burst duration, apparent energy release and 
spectra. 

If our model is a correct representation for short-GRBs 
engines, most of these engines would reside in globu- 
lar cluster s (GCs) where many LMXBs are seemi ngly 
found (Bogdanov et al.|[200l iCamilo fc Rasiol 12001 . It 
also implies that short GRBs should be associated with 
early-type galaxies (no spiral arms) known to form 
their LMXBs in GCs. There is some evidence that 
shor t GRBs reside in th e outskirts of early-type galax- 
ies (jFoxfcRominel [20071) where many GCs are located. 
Finally, the WD ablation in our model provides a po- 
tential unifying framework for the central engines of 
type I Supernovae. In particular, systems observed at 
high latitudes (i.e. > H/R) when the QN goes off, 
could find some interesting applications in the context 
of unu sual type la SNe (e.g. SN2005E; iPerets et al] 
20101)1 and/or sublun i inous type la SNe in general 



Gonzalez-Gaitan et al] (|2010D ). 



7. CONCLUSIONS 
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